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Abstract. We study a family of equations defined on the space of tensor densities 
of weight A on the circle and introduce two integrable PDE. One of the equations 
turns out to be closely related to the inviscid Burgers equation while the other has 
not been identified in any form before. We present their Lax pair formulations and 
describe their bihamiltonian structures. We prove local wellposedness of the corre- 
sponding Cauchy problem and include results on blow-up as well as global existence 
, | of solutions. Moreover, we construct "peakon" and "multi-peakon" solutions for 

d all A 7^ 0, 1, and "shock-peakons" for A = 3. We argue that there is a natural 

geometric framework for these equations that includes other well-known integrable 
equations and which is based on V. Arnold's approach to Euler equations on Lie 
groups. 
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Appendix C. Conservation of H2 for /iDP 
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1. Introduction 



Integrability of an infinite-dimensional dynamical system typically manifests itself 
in several different ways such as the existence of a Lax pair formulation or a bihamil- 
tonian structure, the presence of an infinite family of conserved quantities or at least 
the ability to write down explicitly some of its solutions. In this paper we introduce 
and study two nonlinear partial differential equations and show that they possess all 
the hallmarks of integrability mentioned above. The first of these equations we shall 
refer to as the [iBurgers ([iB) equatio'^ 

(/xB) fi(u t ) - u txx - 3u x u xx - uu xxx = 0, 

and the second as the fiDP equation 

(//DP) ixiut) - u txx + 3/j,(u)u x - 3u x u xx - uu xxx = 0, where = / udx. 

Jo 

Here u(t, x) is a spatially periodic real-valued function of a time variable t and a space 
variable x £ S 1 ~ [0, 1). Both of these equations belong to a larger family that also 
includes the Camassa-Holm equation |CH] (see also |FF] ) 

(CH) u t — u txx + 3uu x — 2u x u xx — uu xxx = 0, 

the Hunter-Saxton equation |HuS] 

(HS) utxx + 2u x u xx + uu xxx = 0, 

the Degasperis-Procesi equation [DPJ 

(DP) u t - utxx + 4uu x - 3u x u xx - uu xxx = 0, 

as well as the /^-equation which was derived recently in jKLMj (we will refer to it as 
the ;uCH equation in this paper 

(//CH) 2fi(u)u x + 2u x u xx + uu xxx = 0, 

all of which are known to be integrable. 

One of the distinguishing features of the CH and DP equations which makes them 
attractive among integrable equations is the existence of so-called "peakon" solutions. 
In fact, CH and DP are the cases A = 2 and A = 3, respectively, of the following family 
of equations 

(1.1) u t - u txx + (A + l)uu x = \u x u xx + uu xxx , A G Z, 

with each equation in the family admitting peakons (see |DHH| ) although only A = 2 
and A = 3 are believed to be integrable (see |DP] ) . One of our results will show that 



each equation in the corresponding fi- version of the family (1.1 ) given by 

(1.2) fi(u t ) - u txx + \/J.(u)u x = Xu x u xx + uu xxx , AgZ, 

also admits peakon solutions. The choices A = 2 and A = 3 yield the /iCH and //DP 



equations, respectively, and as with ( 1.1 ), we expect that these are the only integrable 



^This equation is mentioned in Remark 3.9 of [HoS| and Remark 3.2 of [Luj as the high-frequency 



limit of the Degasperis-Procesi equation, see ( DP I below. Our terminology will be explained in 
Sectio n [Tj 

2 In [KLM] the authors refer to /iCH as the /iHS equation. 
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members of the family (1.2). Moreover, we will show that the /iDP equation admits 



shock-peakon solutions of a form similar to those known for DP, see |Lu . 

In Section [2j we present a natural setting in which all the equations above can be 
formally described as evolution equations on the space of tensor densities (of different 
weight A) over the Lie algebra of smooth vector fields on the circle. In Section [3] 
we present Lax pairs for /iDP and /iB, establishing their integrability. In Section 
[4] we describe the Hamiltonian structure of the equations in (1.2). In particular, we 
consider the bihamiltonian structure of the /iDP equation together with the associated 
infinite hierarchy of conservation laws. In Section [5] we study the periodic Cauchy 
problem of /iDP; we prove local wellposedness in Sobolev spaces and show that while 
classical solutions of /xDP break down for certain intial data, the equation admits 
global solutions for other data. In Section [6] we construct multi-peakon solutions of 
(1.2) as well as shock-peakons of /iDP. In Section [7] we discuss the /iB equation and 
its properties and present a geometric construction to explain its close relation to the 
(inviscid) Burgers equation. Finally, in Section [8] we consider generalizations of the 
/iCH and /iDP equations to a multidimensional setting. 

Left open is the question of physical significance, if any, of the two equations //DP 
and /iB. It is possible that they may play a role in the mathematical theory of water 
waves (as CH, see |CHj ) or find applications in the study of more complex equations 
(as e.g. HS, see |ST] ). We do not pursue these issues here. 

Our approach draws heavily on [KM] and [KLM j and the present work is in some 
sense a continuation of those two papers. 



2. A FAMILY OF EQUATIONS ON THE SPACE OF A DENSITIES 

Perhaps the simplest way to introduce the equations that are the main object of 
our investigation is by analogy with the known cases. We shall therefore first briefly 
review Arnold's approach to the Euler equations on Lie groups and then describe a 
more general set-up intended to capture the /iDP and the //B equations. 

In a pioneering paper Arnold (A] presented a general framework within which it 
is possible to employ geometric and Lie theoretic techniques to study a variety of 
equations (ODE as well as PDE) of interest in mathematical physics (see also more 
recent expositions in |KMj or |KWj ) . Arnold's principal examples were the equations 
of motion of a rigid body in M 3 and the equations of ideal hydrodynamics. 

The formal set-up is the following. Consider a possibly infinite-dimensional Lie 
group G (the "configuration space" of a physical system) with Lie algebra q. Choose 
an inner product (•, •) on q (essentially, the "kinetic energy" of the system) and, using 
right- (or left-) translations, endow G with the associated right- (or left-) invariant 
Riemannian metric. The motions of the system can now be studied either through the 
geodesic equation defined by the metric on G (equivalently, the geodesic flow on its 
tangent or cotangent bundles) or else directly on the Lie algebra q using Hamiltonian 
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reduction^ The equation that one obtains by this procedure on g is called the Euler 
(or Euler- Arnold") equation. Using the inner product it can be reformulated as an 
equation on the dual algebra g* as follows. 

Let (•,•) be the natural pairing between g and g* and let A : g — > g* denote 
the associated inertia operatoiQ determined by the formula (Au, v) = (u, v) for any 
u,v £ g. The Euler equation on g* reads 

(E) mt = — ad^_i m m, m = Au £ g*, 

where ad* : g — > End(g*) is the coadjoint representation of g given by 

(2.1) (ad*m, v) = — (m, [u,v]) 

for any u, v £ g and m £ q* . ad* is the infinitesimal version of the coadjoint action 
Ad* : G x g* — ► g* of the group G on the dual algebra g*. 

In our case G will be the group Diff (S 1 ) of orientation-preserving diffeomorphisms 
of the circle whose Lie algebra g is the space of smooth vector fields T e Diff(5 1 ) = 
vect(S' 1 ). The dual vect*(S' 1 ) is the space of distributions on S* 1 but we shall consider 
only its "regular part" which can be identified with the space of quadratic differentials 
3~2 = {m(x)dx 2 : m S C°°(S' 1 )} with the pairing given by 

[mdx 2 ,vd x ) = / m(x)v(x)dx 
Jo 

(see e.g. Kirillov |KJ). The coadjoint representation of vect(S' 1 ) on the regular part of 
its dual space is in this case precisely the action of vect(5 1 ) on the space of quadratic 
differentials. By a direct calculationm using (2.1) we have 

(2.2) ad* u Q x mdx 2 = (um x + 2u x m) dx 2 
and the Euler equation (|E]) on g* takes the form 

(2.3) m t = — ad*^_i m m = —um x — 2u x m, m = Au. 



Equivalently, when rewritten on vect(S' 1 ), equation (2.3) becomes 
(2.4) Au t + 2u x Au + uAu x = 

so that with an appropriate choice of the inertia operator A (which is equivalent to 
picking an inner product on T e Diff) 



(2.5) 



1 - d 2 for CH, 
fi — d 2 for ^CH, 
-dl for HS, 



one recovers the CH, //CH and HS equ ation^] as Euler equations on vect(S' 1 ); these 
results can be found in [Ml], |KM] and |KLM~} . 

We now want to extend this formalism to include the DP equation as well as 
the two equations of principal interest in this paper: ^iDP and ^B. Admittedly, our 
construction does not have the same beautiful geometric interpretation as that of 
Arnold's and hence perhaps is not completely satisfactory. 



We will be making use of both structures in this paper. 

In infinite-dimensional examples A is often some positive-definite self-adjoint pseudodifferential 
operator acting on the space of smooth tensor fields on a compact manifold. 

^Note that in the case of the HS equation the inertia operator is degenerate, see KM , Section 4, 
for details. 
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The main point is to replace the space of quadratic differentials with the space of 
all tensor densities on the circle of weight Ajj Recall that a tensor density of weight 
A > (respectively A < 0) on S 1 is a section of the bundle (£) X T*S 1 (respectively 
®~ X TS 1 ) and set 

J-\ = \m(x)dx x : m{x) G C 00 ^ 1 )} . 

There is a well-defined action of the diffeomorphism group Diff (S 1 ) on each density 
module T\ given by 

(2.6) T\ 3 mdx x -^mo£ (d x £) X dx X G J^ A , 

which naturally generalizes the coadjoint action Ad* 

space of quadratic differentials. The infinitesimal generator of the action in (|2.6|) is 
easily calculated, 



f G Diff(5 1 ) 
Diff(5 1 ) -> Aut(j^ 2 ) on the 



(2.7) 



L X Q x (mdx x ) = (um x + \u x m) dx x , 



and can be thought of as the Lie derivative of tensor densities. It represents the action 
of vect(S' 1 ) on T\ which for A = 2 coincides with the (algebra) coadjoint action on 
T 2 (that is L 2 uQ 



If we think of (2.7) as defining a vector field on the space T\ then we can consider 



the equation for its flow 
(2.8) 



-um x — Xu x m 



in analogy with (2.3). The substitution m = Au transforms (2.8) into an equation 



on the space of quadratic differentials 



(2.9) 



Aut + Xu x Au + uAu x = 



which is the A- version of (2.4). 

Setting A = 3 and choosing suitable inertia operators A as above we obtain the 
DP, //DP and fiB equations. More precisely, the specific choices of A are in parallel 
with those made for the A 



(2.10) 



2 family in (2.5), that is 

-dl for DP, 
! - d 2 x for //DP, 
-d^ for /jB. 



A 



(1.2). 



More generally, letting A = \i — in (2.9) for any A G Z, we find the equations in 



3. Lax pairs 

An elegant manifestation of complete integrability of an infinite dimensional dy- 
namical system is the existence of a Lax pair formalism. It is often used as a tool 
for constructing infinite families of conserved quantities. This formalism has been 
quite extensively developed in recent years for the CH and DP equations, see e.g. 
|BSS| ICM1 |DHH[ ILS] . Our next result describes the Lax pair formulations for the 
/iDP and //Burgers equations. 



Theorem 3.1. The \xDP and the \iB equations admit the Lax pair formulations 

l^xxx = -Xmip, 

\lpt = -\i>xx - Ulp x + U x ll), 



(3.1) 



^We refer to [O], |OT) or |GR| for basic facts about the space of tensor densities. 
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where A G C is a spectral parameter, ip(t, x) is a scalar eigenfunction and m = Au 
with A = [i — d 2 or A= —d 2 as defined in (2.10). 



Proof. This is a straightforward computation which shows that the condition of com- 
patibility {$t)xxx = (' l Pxxx)t of the linear system (3.1) is equivalent to the equation 
/iDP or /iB when m is given by m = n(u) — u xx or m 



respectively. 



□ 



4. HAMILTONIAN STRUCTURES AND CONSERVED QUANTITIES 

In this section we describe the bihamiltonian structure of //DP (see Section [7] for 
the bihamiltonian structure of the //Burgers equation) as well as the Hamiltonian 



structure for the family of //-equations (1.2) for any A G Z. In the last subsection we 
consider the orbits of Diff (S 1 ) in the space of tensor densities of weight A and give a 
geometric interpretation of one of the conserved quantities. 

4.1. Bihamiltonian structure of //DP. Recall that DP admits the bihamiltonian 
formulation |DHH| 

T 5H j 5H 2 

m = Jo- — = — , 
dm dm 

where 

9 r . . i 



Ho = — / mdx and Ho = — u dx, 
2J 6J 

and the Hamiltonian operators are 

J = m^dxm 1 / 3 (d x - dl)~ l m 1 / 2 d x m 2 / z and J 2 = d x (4 - Sg) (l - d 2 x ) 

with m = (1 — d 2 )u. 

Similarly, the //DP equation admits the bihamiltonian formulation 

_ 5H 5H 2 
m t = Jo^ — = — , 
dm dm 

where now the Hamiltonian functionals Hq and H 2 are 

(4.1) H = J mdx and H 2 = - J ^{u){A^d x uf + ^u 3 ^ dx, 

the operators Jo and J 2 are given by 

J = -m 2 / 3 d x m 1/3 d~ 3 m 1 ^d x m 2/3 and J 2 = -d 3 A = d% 

and m = Au with A = /i — d 2 . The fact that Jo and J 2 form a compatible bihamil- 
tonian pair is a consequence of Theorem 2 in [HW_. It can be verified directly that 
Hq and H 2 , as well as H\ = g / u 2 dx are conserved in time whenever u is a solution 
of //DP (see Appendix [C] for details of this calculation in the case of H 2 ). Using the 
standard techniques we can now construct an infinite sequence of conservation laws 

. . . H-i, Hq, Hi, H 2 , . . . 

see Figure [2j As in the case of the CH and DP equations the above conservation 
laws H n are nonlocal and not easy to write down explicitly for n > 3, while they are 
readily computable recursively in terms of m and its derivatives for negative n. In 
fact, the first negative flow in the /xDP hierarchy is 

/xDP_i :m t = - 2 ( 6160m^ - 13200mm X00 ml + 3600m 2 mlm xxx 



729m 17 / 3 

675m 2 (mm xxxx — 8m 2 x } m x + 27m 3 (3mm xxxxx — 50m xx m xxx ) 
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^DP 3 




Figure 2. Recursion scheme for the [iDP equation. 



4.2. Hamiltonian formulation of ^-equations. Although we expect that /iCH 
and ^DP are the only equations among the family of //-equations in ( 1.2 ) which admit 
a bihamiltonian structure we can nevertheless provide one Hamiltonian structure for 
any A ^ 1. Indeed, if we set 



(4.2) 

and 

(4.3) 



Jo 



1 



tjK + \md x )d x 3 ((A - l)m x + Xmd x ) 



A- 1 



m dx 



then equation (1.2) is equivalent to 



SHq 

5m 



In order to see that Jq is a Hamiltonian operator we first rewrite it in the form 



-m 



(A-l 



)/x d x 



m 



/ X dZ z m l ' X d r m^l x 
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and refer to Theorem 1 of |HW] . 

4.3. Orbits in J-\. From Section 2] w e know that for any value of A a solution m of 
the Cauchy problenj^Jfor equation ( 1.2 ) with initi al data rn-o belongs to the orbit of mo 
in T\ under the Diff (5 1 )-action defined in (2.6). It is therefore tempting to exploit 
this setting to study the family in (1.2) as in the case of the coadjoint action when 
A = 2. We present here two results in this direction leaving a detailed investigation 
of the geometry of the corresponding orbits for general A for a future work. 

The first proposition shows the conservation of H—\ = Jq \m\ l ^dx under the flow of 



( 1.2 ). We let L^(mdx x ) = mo£(d x $,) x dx denote the action of an element £ E Diff (S 1 ) 
on mdx x in J-\. 

Proposition 4.1. The map 



m dx 



H-x[m] = [ \m\ 1/x dx : J-\ 
Jo 



is invariant under the action (2.6) o/Diff(S' 1 ) on T\. 



Proof. This is a straightforward change of variables 

[L £ (mdx A )] = f mo^d^y 
Jo 



H- 



dx 



rn 



1/A 



since £ is a smooth orientation-preserving circle diffeomorphism and m dx x E JF\ 



□ 



In many respects the classification of orbits in T\ resembles that of coadjoint orbits 
of Diff(S' 1 ), cf. [GRJ. In order to state the next result we denote by M.\ the set of 
those elements mdx x in T\ such that m(x) ^ for all x £ S 1 . 

Proposition 4.2. The orbit space (S 1 ) is in bijection with the set M\{0}. 

More precisely, the map 

(4.4) mdx^ i — > sgn(m)(i?_i) A : M.\ — > M\{0}, H-i = [ \m\ 1/x dx, 



is constant on the Diff(5 1 ) -orbits in M\ C J-\ and induces a one-to-one correspon- 
dence between orbits in A4\ and M\{0}. 



Proof. That the map is constant on the orbits follows from Proposition 4.1 To see 
that it is surjective we note that for any nonzero real number a the element adx x 
is mapped to a. In order to prove injectivity we observe that the orbit through an 
arbitrary element mdx x of A4\ contains an element of the form adx x where a is a 
constant. Indeed, given mdx x we define £ E Diff (S 1 ) by 

*(*) ' 



H- 



\m\ l > x dx. 



i Jo 



Then 



L £ (sgn(m)^ lf ix A ) = sgn(m)H x 1 (d x O X dx > 



mdx 



Thus, the action of £ 1 maps m to sgn(m)i? A 1 da; A where sgii(m)H x 1 £ ^\{0} is 
constant. □ 

The stabilizer in Diff(5 1 ) of a constant element adx x in M.\ (a G K) consists 
exactly of the rigid rotations and is isomorphic to S 1 . It follows from Proposition |4.2| 
that the stabilizer of any point mdx x E M.\ is conjugate to S 1 within Diff(5 1 ). In 
particular, each orbit in Ai\ is of the form Diff(5 1 )/5 1 . 



See Section 5.4 for results on wellposedness of the equations ( |1.2[ ) for any A. 
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5. The Cauchy problem for the //DP equation 

In this section we turn to the periodic Cauchy problem for //DP. Rather than 
aiming at the strongest possible theorems we present here only basic results that 
display the interesting behavior of solutions. 

We start with local existence, uniqueness and persistence theorems for the //DP 
equation in Sobolev spaces. Next, we describe the breakdown of smooth solutions 
and prove a global existence result for the class of initial data with non-negative 
momentum density. At the end of this section we briefly discuss local wellposedness 
of the whole family ( |1.2[ ) of /i-equations. 

In order to prove these results we will need Sobolev completions of the group of 
circle diffeomorphismg^j Diff (S 1 ) and its Lie algebra of smooth vector fields vect(S' 1 ). 
We denote by H s = H S (S 1 ) the Sobolev space of periodic functions 

H^S 1 ) = {/; = ^2v(n)e 2winx : \\vf H . = ^ \A^v(n)\ 2 < 00} 

n n 

where s 6 M and the pseudodifferential operator A s = (1 — d 2 ) 3 ^ 2 is defined by 

A*v(n) = (1 + 4vrV) s/2 /)(ra). 
In what follows we will also make use of another elliptic operator 

(5.1) A 2 : H'(S X ) - H'-^S 1 ), Ajt, = //(<;) - v xx 
whose inverse can be easily checked to be 

(5.2) A~ 2 v(x) = ( y - I + ^ J v(x) dx+(x-^J J v{y) dydx 



v(z) dzdy + / / / v(z) dzdydx. 
Jo Jo Jo Jo 

5.1. Local wellposedess in Sobolev spaces. We begin with the //DP equation. 
Our proof follows with minor changes the approach of |M2 to the CH equation based 
on the original methods in [EM] developed for the Euler equations of hydrodynamics. 
We include it for the sake of completeness. The two above papers and [KLMJ will be 
our main references for most of the basic facts we use in this section. 
We can write the Cauchy problem for //DP in the form 

(5.3) u t + uu x + 3//(m) d x A~ 2 u = 

(5.4) u(0) = n . 



Note that applying A 2 to both sides of (5.3) we obtain the equation in its original 
form given in the Introduction. 

Theorem 5.1 (Local wellposedness and persistence). Assume s > 3/2. Then for 
any uq £ H S (T) there exists a T > and a unique solution 

u e C({-T,T),H S ) nC 1 ((-T,T),F s - 1 ) 



of the Cauchy problem ( |5.3| )-( |5.4[ ) which depends continuously on the initial data uq. 
Furthermore, the solution persists as long as Ollc 1 stays bounded. 



^Although no longer a Lie group, Diff s (S 1 ) retains the structure of a topological group for a 
sufficiently high Sobolev index; see below. 
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Our strategy will be to reformulate (|5.3[ )-(5.4 ) as an initial value problem on the 
space of circle diffeomorphisms DifP(S' 1 ) of Sobolev class H s . It is well known that 
whenever s > 3/2 this space is a smooth Hilbert manifold and a topological group. 
We will then show that the reformulated problem can be solved on Diff s (S' 1 ) by 
standard ODE techniques. 

Let u = u(t, x) be a solution of ^DP with initial data uq. Then its associated flow 
t — > x), i.e. the solution of the initial value problerrjj 

(5.5) £(t,x) = u(t,£(t,x)) t (,(0,x)=x 

is (at least for a short time) a C 1 smooth curve in the space of diffeomorphisms 
starting from the identity e G Diff s (S' 1 ). Differentiating both sides of this equation 



in t and using (5.3) we obtain the following initial value problem 



(5.6) 1 = -3 M (£ o C 1 ) {d x A~\i o r 1 ) J o £ =: -F(£, £) 

(5.7) Z(0,x) = x, i(0,x) = u (x) 

In this section we make repeated use of a technical result which we state here for 
convenience (for a proof refer to |M2 , Appendix 1). 

Lemma 5.2. For s > 3/2 the composition map £ — » cj o £ mitt an -fP function u 
and the inversion map £ — > £ _1 are continuous from Diff s (S' 1 ) to H S (S 1 ) and from 
Diff s (S' 1 ) to itself respectively. Moreover, 

(5.8) ||u,o£||tf S <<7(l + ||£||^)H|^ 

with C depending only on inf \d x £\ and sup \d x £\- 

It will also be convenient to introduce the notation = o A o R^-i for the 
conjugation of an operator A on ^(S" 1 ) by a diffeomorphism £ £ Diff s (5 1 ). The 



right hand side of (|5.6|) then becomes 

— /. 



(5-9) F(Z,£)=A-]d x ,th(U) 



where h(£, uj) = 3a; J oj o £ 



Proof of Theorem \5.1\ From Lemma 5.2 we readily see that F maps into H 3 ^ 1 



We aim to prove that F is Frechet differentiable in a neighborhood of the point 
(e,0) in DifF x H s (S l ). To this end we compute the directional derivatives d^Fr^) 
and duFt£ iU) ) an d show that they are bounded linear operators on H s which depend 
continuously on £ and u. We use the formulas 

(5-10) ^A-|( V ) = -A-| hf^^A^A^ 

(5.H) 9fy«) = [« o r 1 9»A] £ 

and 

(5.12) d(:h(£ >U j-) (v) = 3 u / a; o £ i ~ 1 d x (v o £ _1 ) da; 



Here "dot" indicates differentiation in t variable. 
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along with (5.9) to get 

d 6 F (M (v) = 3 { - v o e _1 A; 2 ^(u; o C 1 ) f u> o C^x 

+ A- 2 a x ((v o r 1 )^^ o r 1 )) y (w o r 1 ) 
- A; 2 a x (o; o r 1 ) y (w o r 1 ) ^,(« o r 1 ) ^} ° e 

which with the help of the identity 

(5.13) A; 2 ^ = -l + jU 

can be simplified to 

d^F^^(v) = 3 <.vu J uj o £ _1 <ix ~ v (^J u C^dx) 

+ A~^d x ^(v d x ^u) J ujo£- l dx 

- A~£d X £U J u o C _1 5 2 (?; o f -1 ) da; j . 

On the other hand we have 

d LO h^ ^(v) = 3v / u) o dx + 3uj / u o £~ da; 



and similarly 

d w Ffa u )(v) = —3A~£d x £V J uj o £~ l dx - 3A~^d x ^u j vo^dx. 

In order to show that v — * d^F^^(v) is a bounded operator on H s it is sufficient 
to estimate the sum 

\\d^ >uj) (v)\\ L2 + Wd^F^iv^oC^Hs-i- 



Using Cauchy-Schwarz, Lemma |5.2| and the formulas above both of these terms can 
be bounded by a multiple of IMI 2 ^ which gives the estimate 

ll<%%^)0)||^ < C?||w|| jE r- ||^||lr« 

where C depends only on inf \d x £\ and sup \d x £\. A similar argument yields the bound 
for the other directional derivative 

\\duF(£ tC0 ){v)\\ H s < C|MliHMI# a - 

It follows that F is Gateaux differentiable near (e, 0) and therefore we only need to 
establish continuity of the directional derivatives. Continuity in the w-variable follows 
from the fact that the dependence of both partials on this variable is polynomial. It 
thus only remains to show that the norm of the difference 



is arbitrarily small whenever £ is close to the identity e G Diff s (S' 1 ), uniformly in v 
and oj G H S (S 1 ). Once again, it suffices to estimate the terms 

(5.14) ptF^W-dtF^WWv + Wd^F^v) -a^)(«))L-i 
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We proceed as above. Using formula (5.2) and Lemma 5.2 together with the fact 
that for any r > we have a continuous embedding H r (S 1 ) ^ C r ~ 1 / 2 (S 1 ) with a 
pointwise estimate 



\\g(x) - g(y)\\ < WgWuAx-y 



ir— 1/2 



we can bound both terms in ( 5.14| ) by 



C\\co 



|2 



\ v \\H° 



+ w 



-3/2 



duF(p^\ follows from a similar estimate. 



Continuity of £ 

In this way we obtain that F is continuously differentiate near (e, 0) and applying 
the fundamental ODE theorem for Banach spaces (see e.g. Lang (Lj) establish local 
existence, uniqueness and smooth dependence on the data kq of solutions £(t) and 



£(t) of (5.6)-(5.7). Local wellposedness of the original Cauchy problem ( |5.3| )-(5.4) 
follows now at once from u = £, ° an d the f a °t that DifP is a topological group 
whenever s > 3/2. 

In order to complete the proof of Theorem 5.1 we need to derive a C 1 bound on 
the solution u. The standard trick is to use Friedrichs' mollifiers J e with (0 < e < 1) 
to derive the inequality 
-l 



d 
dt 



\JfU\ 



[ A s J e ud x A s J e v 2 dx-^fi(u) [ A 2 J € ud x A s A 2 J e udx 
Jo 2 J 



< 



\u\\ci\\u\\'hs + ||^a;A S A 



-2 



J £ u\\ L2 \\A s J € u\\ L2 < 



\u\ 



c 



1 \\u 



[2 



where the first estimate is standard (see e.g. jT], Chapter 5) and the second follows 
easily with the help of (5.1) and (5.13). Passing to the limit with e — * and using 
Gronwall's inequality yields the persistence result. □ 



5.2. Blow-up of smooth solutions. A family of examples displaying a simple 
finite-time breakdown mechanism of a /^DP solution is described by the following 
result. 

Theorem 5.3. Given any smooth periodic function uq with zero mean there exists 
T c > such that the corresponding solution of the fj,DP equation stays bounded for 
t <T C and satisfies \\u x (t)\\oo /* °o as t f T c . 

Proof. Let u be the solution of //DP with initial data uq. Differentiating the equation 
with respect to the space variable we obtain 

utx + uuxx + u 2 x + 3[i(u)d 2 A~ 2 u = 0. 

Since d 2 A~ 2 = [i — 1, using conservation of the mean, we find that 

(5.15) u tx + uu xx + u 2 = 3fi(u) (u - n(u)) = 



by the assumption on the data Uq. If we let £(t) denote the flow of u as in (5.5) then 

d x i = u x °i dxi 

and setting w = d x £/d x £ we find that 

dxjdxt-m) 2 f ^ w 

W t = /a = \ u xt + UxxU) o f. 

With the help of these formulas we rewrite ( |5.15[ ) in the form 

wt + w =0 
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and, since w(0,x) = uq x (x), solve for w to get 

_ 1 

t+(l/uo x (x)) 

Furthermore, our assumptions are such that it is always possible to find a point 
x* G S 1 such that uq x (x*) is negative. Setting T c = —1/uq x {x*) we conclude that the 
solution must blow-up in the L°° norm as expected. □ 

More sophisticated break-down mechanisms for //DP can be demonstrated but we 
will not pursue them in this paper. 

5.3. Global solutions. Our global result for //DP is obtained under a sign assump- 
tion on the initial data of the type that has been used previously in the studies of 
equations such as CH, /iCH and DP. What follows is again a basic result. 

Theorem 5.4. Let s > 3. Assume that uq G H S (S 1 ) has non-zero mean and satisfies 
the condition 

Ajuo > (or < 0). 

Then the Cauchy problem for fiDP has a unique global solution u in C(R, H S (S 1 )) H 



Proof. From the local wellposedness and persistence result in Theorem 5.1 we know 
that the solution u is defined up to some time T > and that in order to extend it 
we only need to show that the norm ||u x (i, -)||oo remains bounded. 

On the one hand, given any periodic function w(x) differentiating the formula in 



(5.2) we readily obtain the estimate 



(5-16) \\d x w\\oo< ||A>|| L1 . 

On the other hand, recall from Section [2] that for any solution u of //DP with 
initial condition uq G H s (S 1 ) the expression m = A^u is a tensor density in J-^S ) 



and hence from the Diff s (S' 1 )-action formula in (2.6) we get the following (pointwise) 
conservation law 

(5.17) — ({A^u) o f (d x tj) 3 ^ = ( - u txx - uu xxx + Sfx(u)u x - 3u x u xx ) o £ (d x £) 3 = 0, 

where t — ► £(£) G Diff s (S' 1 ) is the associated flow (i.e. £ = u o £ and £(0) = e). Thus, 
for any x G S 1 as long as the solution exists it must satisfy 

Aju(i, C(t, x)) (a x e(t, x)f = Alu (x) 
and therefore by the assumption on no it follows that 

A 2 ^u(t, x) > for any x G S 1 
as long as it is defined. The above inequality together with the conservation of the 



mean and (5.16) gives 

«l 



\\d x u(t, ^Woo <\\A 2 u(t)\\ L1 = / A 2 u(t, x) dx = u(u ) < oo 

Jo 

and we conclude that the solution u must persist indefinitely. □ 
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5.4. Local wellposedness for any A. For completeness we state here a local well- 



posedness result for the family (1.2) of ^-equations for any A. Using the substitution 
v = Av = A^u and the fact that, as in the case of /xDP, the mean of the solution u is 



conserved, we can write the Cauchy problem for (1.2) in the form 



(5.18) 



u t + uu x + A d x [ Xfi(u)u + 



A 



-u, 



(5.19) 



u(0) = u . 



Theorem 5.5 (Local wellposedness). Let s > 3/2. For any uq 6 H S (S 1 ) there exist 
a T > and a unique solution 



£ C((-T,T),H s )nC 1 ((-T,T),H 



s-l\ 



of (5.18)-(5.19) which depends continuously on the initial data uq. Furthermore, the 



solution persists as long as \\u(t, stays bounded. 



The proof of this theorem is analogous to Theorem 5.1 and we omit it. 
In the previous section we proved a global existence result for /iDP. Here, we 
formulate a similar result for any positive A. 



For a fixed A > let u = u(t,x) be any smooth solution of the problem (5.18)- 



(5.19) and let £(i) £ Diff s (S' 1 ) denote the corresponding flow of diffeomorphisms of 
the circle. As before, the argument rests on the pointwise conservation law obeyed 
by the solutions. Namely, using the formula for the action of Diff s (S' 1 ) on A-densities 



T\(S ) in (2.6) we find the analogue of (5.17) to be 



d_ 
dt 



(A^u) o £ (d x C) x ^j = (A^u t - uu xxx + \^{u)u x - \u x u xx ) o £ (d x £,)' 



so that as long as u(t, x) exists 

Alu(t,£(t,x))(d x at,x)) X = Aluo(x), xeS 1 
and therefore we have 

Theorem 5.6. Fix A > 0. Assume that u$ S H s (S l ) with s > 3 has non-zero mean 
and satisfies the condition 

A^uq > (or < 0). 



Then the Cauchy problem (5. 18) -(5. 19) has a unique global solution in C (M, H s (S 1 )')C\ 
C^R,^- 1 ^ 1 )). 



6. Peakons 

The CH and DP equations famously exhibit peakon solutions, see e.g. |CH|IDHH] . 
In this section we show that each of the /^-equations ( 1.2 ) parametrized by A S Z also 
admits peakons. Let us point out here that the fact that /uCH admits peaked solutions 
went unnoticed in [KLMJ, so that the present discussion provides an extension of the 
results of that paper. 

We will investigate existence of traveling- wave solutions of the form u(t, x) = f(x — 
ct). This will lead us to expressions for the one-peakon solutions. Subsequently, we 
will analyze the more general case of multi-peakons. 
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6.1. Traveling waves. We first consider equation (1.2) where fi(u) is replaced by a 
parameter v 6 R, i.e. 

(6.1) - u txx + \vu x = \u x u xx + uu xxx . 

After analyzing traveling-wave solutions of this equation we will impose the two con- 
ditions 

jx{u) = v and period(ii) = 1 



to find which of these are traveling waves of equation ( 1.2 ). We assume that A 7^ 0, 1. 
Substituting u(t, x) = <p(x — ct) into equation (6.1 ) we find 

cip xxx + \v<fx = Xip x ip xx + W xxx ■ 
We write this in terms of / := ip — c and integrate to find 



(6.2) 



Xvf + a 



A-l 



fi + ffz 



where a in an integration constant. We multiply equation (6.2) by 2f x 2 f x and 

2vf x + 



integrate the resulting equation with respect to x. The outcome is 

2a 



x _ 1 f X ~ 1 + b = f A - 1 f x , 
where b is another constant of integration. This equation can be written as 



A-l f 2 



(6.3) 



2v{if> - c) + 



2a 
A-l 



+ 



A solution of this equation yields the shape of a traveling wave within an interval 
where ip is smooth, and gluing such smooth wave segments together at points where 
ip = c produces the complete collection of traveling wave solutions of (6.1). The 
analysis proceeds along the same lines and yields similar results as for CH (or as for 
HS if v = 0) cf. |Le| . Here we choose to consider the peaked solutions. 

For completeness we point out that the equation analogous to (6.3) when A = or 
A = 1 is 



(6.4) 



-2a + b(<p-c), A = 0, 

2v(ip - c) + 2aln|v? - c| + b, A = 1. 



6.2. Periodic one-peakons. The periodic peakons arise when b = in (^Pl 
Indeed, setting 6 = and replacing a by a new parameter m via a = (A — l)(c — m)v 
in (6.3), we find 

<P 2 X = 2v{ip-m). 
Solving this equation, we find the peakon solution u(t, x) 



(p(x — ct) where 



(6.5) 



(p(x) 



m + —x 



for x € 



2(c 



m) 



2(c 



mj 



and ip is extended periodically to the real axis. The parameters are required to satisfy 
(c — m)/v > 0. The period of 99 is fixed by the condition that ip = c at the peak. 

Equation (6.5 ) defines a periodic peaked solution of (6.1 ) for any choices of the real 
parameters m, c, v. In order to determine which of these solutions are solutions with 
period one of the /i-equation (1.2) we have to impose the conditions mean(^) = v 
and period(<£>) = 1. A computation shows that mean(c^) = m + and period(c^) = 



10 



We henceforth assume v 7^ 0. 



16 



JONATAN LENELLS, GERARD MISIOLEK, AND FERIDE TIGLAY 



2(c-m) 



Therefore, u(t,x) = cp(x — ct) where ip is given by (6.5) is a period-one 



solution of equation (j 1 . 2 p if and only if 

c 



2(c — m) 



1. 



m c 
— > 0, and v = m H 



m 



Solving these equations we find m = 23c/26 and v = 12c/13. This leads to the 
following result. 



Theorem 6.1. For any c£M and A ^ 0, 1, equation (1.2) admits the peaked period- 
one traveling-wave solution u(t, x) = ip(x — ct) where 

c 



(6.6) 



ip(x) 



26 



[12x z + 23) 



for if [— g, |] arid * s extended periodically to the real line. 



We note that the one-peakon solutions of (1.2) are the same for any A, that they 
travel with a speed equal to their height, and that fj,(<p) = 12c/13. 

6.3. Green's function. For the construction of multi-peakons it is convenient to 
rewrite the inverse of the operator A 2 = \i — d^. in terms of a Green's function (an 
explicit formula for the inverse A~ 2 was previously given in equation (5.2)). We find 

(A~ 2 m)(x) = / g(x — x')m(x')dx', 



where the Green's function g{x) is given by 
(6.7) g{x) 
and is extended periodically to the real line. In other words, 

x,x' e [0,1) ~ S 1 . 



I 23 

-x(x-l) + — for x£[0,l)-S\ 



, (x — x') 2 \x — x'\ 13 

g(x - x') = H , 

yv ' 2 2 12' 



Note that g(x) has the shape of a one-peakon (6.6) with c = 13/12 and peak located 
at x = 0. In particular, fi(g) = 1. 

6.4. Multi-peakons. The construction of multi-peakons for the family of ^-equations 
( 1.2 ) is similar to the corresponding construction for the family ( 1.1 ) cf. [DHH2]. Just 
like for CH and DP, the momentum m = A 2 u of an iV-peakon solution is of the form 

N 

(6.8) m(t,x) = Y,Pi(tnx-Q i (t)), 

i=l 

where the variables {q l ,Pi}i evolve according to a finite-dimensional Hamiltonian 
system. Indeed, assuming that m is of the form (6.8), the corresponding u is given 
by 

N 



i=l 



We assign the value zero to the otherwise undetermined derivative </(0), so that 



(6.9) 



g'{x) :-- 



0, x = 0, 
x—\, < x < 1. 
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Figure 3. The periodic Green's functions g(x) and g^{x) corresponding to the op- 
erators A = 1 — d x and A = // — d x , respectively, and their difference g(x) — g^x). 



This definition provides a naive way to give meaning to the term —\u x m in the PDE 



(2.8). Substituting the multi-peakon Ansatz (6.8) into the //-equation (1.2) then 



N 



shows that pi and q l evolve according to 

JV 

q* = ^Pjgtf - Q j ), Pi = -(A - i) ^PiPjg'tf - g j )- 

3=1 3=1 

A more careful analysis reveals that this computation is in fact legitimate and yields 
the following result. Note that the peakons fall outside the result of Theorem |5.5 
since for a peakon u ^ H S (S 1 ) for s > 3/2. 



Theo rem 6.2. The multi-peakon (6.8) satisfies the [i-equation (1.2) in the weak form 
(5.18) in distributional sense if and only if {q l ,Pi}\ evolve according to 

(6.10) <f = «(<f), 

Pi = -(A - l)pi{u x (q 1 )} 

where {u x (q 1 )} denotes the regularized value of u x at q l defined by 

N 

(6.H) {u x (q l )}:=J2Pj9'(Q l -q J ) 

3 = 1 



and g'{x) is defined by (6.9). 
Proof. See Appendix [A] 



□ 



The formulas for the peakons of the //-equations (1.2) along with the expressions 



for the previously known peakons of the CH and DP family (1.1 ) are summarized in 



Table 6.4 At any particular time, the multi-peakon is a sum of Green's functions for 
the associated operator A, see Figure [3} Since the inertia operator for the //-equations 
(A = fjL — d%) is different from that of the family (1.1 ) (A = 1 — d%), the new class 
of "//-peakons" have a different form with respect to the CH and DP peakons. For 
a given family of equations the Green's function is the same for all A, but the time 
evolution of the positions q % and momenta pi of the peaks of course depends on A and 



is for both families given by the system (6.10). 



6.5. Poisson structure. The equations (6.10) take the Hamiltonian form 

q i = {q\H }, pi = {pi,H }, 
with respect to the Hamiltonian 



IS 
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Equation family 


Green's function (x £ S 1 ~ [0, 1)) 


Multi-peakon u(t, x) 


( 
( 


1.1 
1.2 


) 
) 


n (~\ - COSk(x-l/2) 

9{X) - 2sinh(l/2) 
5 (x) = i(x-i) 2 + i 


E£ift(*Mx -«*(*)) 



Table 1. The periodic multi-peakon solutions of the family (1.1) 
(which includes CH and DP) and of the corresponding //-family (1.2) 
(which includes /xCH and //DP) when A ^ 0, 1. 



and the Poisson structure 



{Pi,Pj} 



A- 1 



A 



PiPjG 



A-l 



PjG'{q l - q j 



1 

A2 



G(q l -q j ), 



where G(x) = f£ g(x')dx' cf. |DHH2j . In order to make sense of G(q l — q 3 ) although 
G is not periodic, we restrict q l ,q J to be local coordinates on S 1 satisfying say < 
q l ,q 3 < 1. Then —1 < q % — q 3 < 1 and G{q l — q 3 ) is well-defined. This Poisson 



structure can be derived from the Poisson structure corresponding to Jo in (4.2) by 
noting that 

{m(x), m(y)} Jo = - [G{x - y)m x {x)m x {y) - \G'(x - y)m x (x)m(y) 
+ XG'(x — y)m(x)m x (y) — X 2 G"(x — y)m(x)m(y)] . 



Indeed, let mo(x) = YliP0i$( x ~ 1o) be a multi-peakon of the form (6.8). Let fa : 
5 1 — > R (resp. fa : S 1 — > M) be a function which takes the value 1 (resp. x) in a small 
neighborhood of x = q$ and the value elsewhere. Then, for mix) = ^2iPiS(x — q l ) 
sufficiently close to mo, we obtain 



fa(x)m(x)dx = pi, / fa(x)m(x)dx 
is 1 Js 1 

We now find {pi,Pj} from the computation 

{Pi,Pj}=l I fa(x)4> j (y){m(x),m(y)}j dxdy 
Js 1 Js 1 

We then compute 
U l Pi,Pj} 



qPi 



PiPj&'tf-qi). 



s 1 Js 1 
A-l 



fa(x)(pj(y){m(x) ,m(y)} Jo dxdy 
PlPj G'(q l - J) + ( > ^l) 2 q %p j G"(q i 



Since {q l Pi,Pj} = q l {pi,Pj}+{q l ,Pj}pi and we know the value of {pi,Pj}, this gives the 
expression for {q l ,Pj}- Finally, the expression for {q l ,q J } is obtained by calculating 



{q l Pi,q J Pj}= / fa{x)fa{y){m(x),m(y)} Jo dxdy 
Js 1 Js 1 

and comparing the result with the expression for 

{q % Pi, q 3 pj} = q l q 3 {pi,Pj\ + q l {pi, q J }pj + q>{q\Pj}Pi + {q\ q 3 }PiPj 
obtained by substituting in the expressions for {pi,Pj}, {Pi,q J }, and {q l ,Pj}- 
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6.6. Multi-peakons for /iCH. For /iCH (A = 2) the equations (6.10) for {q l ,pi} 
are canonically Hamiltonian with respect to the Hamiltonian 

1 N 



and describe geodesic flow on T 



N . 



o N 



(here T N = S 1 x • • • x S 1 denotes the TV-torus) 



with respect to the metric g a b with inverse given by 



In this case the reduction to multi-peakon solutions can be elegantly understood in 
terms of momentum maps using the ideas of |HM] . Indeed, let S = {1, 2, . . . , TV} be 
a set consisting of N elements. We consider (q l ,pi) as an element of T*Emb(5', S 1 ) ~ 



5 N 



. Then the map 



N 



J sing : T*Emb(,S, S 1 ) -» vectC^ 1 )* : (q\ Pi ) » J2 Pi S(x - q 1 ) 



is an equivariant momentum map for the natural action of Diff (S 1 ) on T*Emb(S', S* 1 ). 
Since momentum maps are Poisson, Using maps the solution curves in T*Emb(5, S 1 ) 
which satisfy (6.10) to solution curves in vect(S' 1 )* which satisfy (2.3). Note that 
h = J* in H is the pull-back of the Hamiltonian H on vect(S' 1 )* defined by H = 



| Jq 1 umdx. 



6.6.1. Periodic two-peakon for fj,CH. In the case of two peaks for the //CH equation 
(N = 2, A = 2), the system (6.10) is given explicitly by 

q 1 =pig(0) +P2g(q 1 - q 2 ), 

<f =pig(q 2 -q 1 )+P2g(o), 

Pi = -P2 = -PW2g'(q 1 - q 2 )- 
Letting Q = q 2 — q 1 and P = P2 — Pi and using that 

13 Hi PlP2 (n . 2 



h = ± + 

24 16 2 



'-(Q-\Q\), 



we find the system 



Q 
P 



l 



^(Q 2 - IQI), 



a 



'Q 2 -\Q\ 



where a = 2h — fs^#- This leads to the following equation for Q: 



12 16 



Q{Q l - 101) = (2Q Z + a(Q z - 101)) I Q- Jsgn(g) ) . 



The solution to this ODE can be expressed in terms of inverses of elliptic functions. 
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6.7. Shock-peakons for //DP. In addition to the peakon solutions the DP equation 
admits an even weaker class of solutions with jump discontinuities called "shock- 
peakons", see |Lu| . We will show here that //DP also allows shock-peakon solutions. 
We will seek the solutions of the form 



N 



mi 



i=i 



i(t, x) = Y^ [Pi(t)S(x - + si(t)8'(x - q\t)) 
The corresponding u is 

(6.12) u = ^(p i g{x-q i ) + s i g'(x-q i ) 



N 



t=l 



Substituting this into the equation mt = —um x — Xu x m, we find formally 

q { = u{q l ), 

Pi = (A - 1) (siU xx (q l ) - PiU x (q 1 )) , 
Si = -(A - 2)siU x (q l ). 

The following theorem states that in the case of A = 3 (this is the case of /uDP), this 
formal computation can be rigorously justified provided that the u x and u xx terms 
are appropriately regularized. The reason that such weak solutions can be m ade sense 
of in the case of //DP is that the term ^^u x in the weak formulation (5.18 ) is absent 
when A = 3. 



Theorem 6.3. The shock-peakon (6.12) satisfies the weak form (5.3) of fiDP in 
distributional sense if and only if {q t ,Pi, S{}^ evolve according to 



(6.13) 



where 



<? =u(q i ), 

Pi = 2{si{u xx {q % )} - Pi{u x {q 1 )}) , 
Si = -Si{u x (q 1 )}, 



N N 

K(<f )} = ^PiSV - Q j ) + Y. s i> 

3=1 3=1 



N 



{u xx (q 1 )} = ^Pj, 

3=1 



and g'{x) is defined by (6.9). 
Proof. See Appendix |Bj 



□ 



7. The //Burgers equation 

In this section we discuss the fiB equation and its properties. We introduced this 
equation in Section [2] as the A = 3 analogue of the Hunter-Saxton equation. Here 
we take a different view and employ Riemannian geometric techniques to display its 
close relationship with the inviscid Burgers equation. 

7.1. Burgers' equation and the L 2 -geometry of DifP — > DiP/S 1 . The con- 
structions in this subsection are well known. First, recall that the inviscid Burgers 
equation 

(B) ut + uu x = 
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is related to the geometry of the weak Riemannian metric on Diff s (S' 1 ) which on the 
tangent space T^DifF s (5 1 ) at a diffeomorphism £ is given by the I? inner product 

(7.1) (V,W) L 2= f V{x)W{x)dx 



o 



where V, W £ TJjDifP. It is not difficult to show that this metric admits a unique 
Levi-Civita connection V whose geodesies rj(t) in DifF(S' 1 ) satisfy the equation^] 

(7.2) Vf,r] = i)= (u t + uu x ) or] = 

and hence correspond to (classical) solutions of the Burgers equation. Here r/(t) is 
simply the flow of u(t, x) so that r){t, x) = u(t, i]{t, x)) and hence the second equality 



in (7.2) follows at once from the chain rule. 

Consider next the homogeneous space DiffQ = DifF / S 1 where S 1 denotes the sub- 
group of rotations. For s > 3/2 this space is also a smooth Hilbert manifold whose 
tangent space at the point [e] in Diff^S 11 ) corresponding to the identity diffeomor- 
phism can be identified with the periodic Sobolev functions with zero mean 



r [e] DifF = H^S 1 ) = Le H^S 1 ) : J* 



w(x) dx = 



It is useful to think of Diff^S" 1 ) as the space of probability densities on the circle and 
view the map ir : Diff s — ► DiffQ as a submersion given by the pull-back [£] = 7r(£) = 
£*(dx) where dx is the (fixed) probability measure with density 1 (see e.g. |K Yj ) . The 
fibre through £ 6 Diff s (S' 1 ) consists of all diffeomorphisms that are obtained from £ 
by a composition on the left with a rotation and thus is just a right coset of the 
rotation subgroup S 1 which acts on Diff s (S' 1 ) by left translations. Furthermore, the 



L metric (7.1) is preserved by this action and tt becomes a Riemannian submersion 
with each tangent space decomposing into a horizontal and a vertical subspace 

T^DifP = ffc(T e Diff a ) L 2 Q f (T 5 Diff s ) 



which are orthogonal with respect to (7.1). The two orthogonal projections P^ 



TgDifP — > T^-^Diffg and : T^DifF — ► IR are given explicitly by the formulas 

(7.3) P^W) = W- f W{x) dx and Qt(W) = [ W{x) dx. 

Jo Jo 



Finally, as with any Riemannian submersioq [ note that a necessary and sufficient 



condition for a curve rj(t) in Diff s (S' 1 ) to be an L 2 geodesic satifying (7.2 ) (and hence 
correspond to a solution of Burgers' equation) is that 

P^V^r) = QnVrfi] = 0. 

We are now ready to introduce the //B equation in this set-up. 

7.2. Lifespan of solutions, Hamiltonian structure and conserved quantities 
of the ^Burgers equation. We first characterise local (in time) smooth solutions. 

Theorem 7.1. A smooth function u = u(t,x) is a solution of the fJ,B equation 

Utxx ~\~ 3lix^xx ~\~ UU X xx — 

if and only if the horizontal component of the acceleration of the associated flow rj(t) 
in Diff s (5' 1 ) is zero i.e. P^V^t) = 0. In fact, given any uo € H S (S 1 ) the flow of u has 
the form rj{t, x) = x + t(uo(x) — Uq(0)) + 7/(t, 0) for all sufficiently small t. 



^In fact, the geodesic equation is readily obtained from the first variation formula for the energy 



functional E(rf) = | /* H'iMll^di of the L 2 metric (O on Diff^S 1 



12 t 



See e.g. O'Neill ON for details on Riemannian submersions. 
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Proof. Integrating the equation with respect to the x-variable we get 

(7.4) u tx + ul + uu xx = \ [u 2 x + uu xx )dx = 

J 

and integrating once again gives 

u t + uu x = c{t) 

where c{t) is a function of the time variable only. Since the mean of a fiB solution 
need not be preserved in time we obtain the equation^ 

(7.5) u t + uu x = n(u t ). 

On the other hand, if rj(t) is the flow of u in Diff s (S' 1 ), that is 

rj(t, x) = u(t, r](t, x)), r](0,x)=x 



then from (I7.2I) and ( 7.3 ) we compute the horizontal component of the acceleration 



of ?7(t) in Diff^S 11 ) to be 

P^nV = V ~ / rjdx = (ut + uu^ or] — I (ut + uux) o r\ dx. 
Jo Jo 

Using these formulas and integrating by parts we find that the equation P^V^?) = 

is equivalent to 

(7.6) (ut + uu x ) o rj = rj = rjdx = rj o r/^ 1 dx 

Jo Jo 

i ,i 

(u t + uu x ) dx = ut dx 
o Jo 

which establishes the first part of the theorem. 

In order to prove the second statement it suffices to observe that from (7.6) we 

have in particular 

rj(t,x)= / f] o r/ _1 (t, x) dx = rj(t, 0) 



which immediately implies 

r](t, x) - rj(t, 0) = x + t(u (x) - u (0)) 
for any < x < 1 and any t for which r\ is defined. □ 



Remark 7.2. It is easy to verify that equation (7.5) is also equivalent to VriP ri r) = 
which can be interpreted as saying that the horizontal component of the velocity of the 
flow rj(t) of the fiB equation is parallel transported along the flow. 

Remark 7.3. A similar construction can be used to derive the (iB equation from the 
right-invariant L 2 metric on Diff ,s (5' 1 ) defined by 

(V, W)(= f V o C\x)W o cHx) dx where V, W G T ? Diff s and f G Diff s . 
Jo 

In this case the corresponding orthogonal projections are P^(W) = W — W d£ = 
W — Q^(W) and the geodesic equation inDiS s (S 1 ) reads V^t) = 7) + 2i]d x r)(d x 'n)~ 1 = 
0. Proceeding as in the proof of Theorem \7.1\ we then get 

= PrjVrjT] = (ut + 3UU X - fJ,(ut)) O 7] 



13 



In the case when the mean n(u) is independent of time (7.5| becomes the standard (inviscid) 



Burgers equation. The presence of the mean is also the reason for our terminology. 
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which after rescaling the dependent variable yields ( 7.5 ) . 

It is therefore not surprising to find that the uB equation shares a number of 
properties with the Burgers equation. 

Corollary 7.4. Suppose that u(t, x) is a smooth solution of the fiB equation and let 
it(0, x) = uq(x). 

(1) The following integrals are conserved by the flow of u 

1 nl 

(u — fJ.(u)) p dx = / (no — fi(uo)) p dx, p = 1,2,3.... 
Jo 

(2) There exists T c > such that ||ux(t)||oo /* °° &s t f T c . 

(3) The \iB equation has a bi-Hamiltonian structure given by the Hamiltonian 
functionals 

Hn = — f m dx, H 2 = — / u 3 dx, 
2J 6J 

and the operators 

J = -m 2 /' i d x m 1 / 3 d~ 3 m 1/3 d x m 2 l 3 , J 2 = d 5 x . 

Proof. The first statement follows by direct calculation. Regarding the second, it 
suffices to consider the equation in (7.4) and apply the argument that was used 
in the proof of Theorem 5.3 The third statement follows from a straightforward 
computation establishing 

7 SH 5H 2 

J - — = J 2 - — = -um x - 3u x m 
dm dm 

where m = —d x u and our earlier observation from section [4] that Jo and J 2 are 
compatible. □ 



8. Multidimensional uCH and uDP equations 

In this final section we briefly consider possible generalizations of the uCH and 
uDP equations to higher dimensions. Let 3t(T n ) be the space of smooth vector fields 
on the n-dimensional torus T n comprising the Lie algebra of Diff(T n ). Let A be a 
self-adjoint positive-definite operator defining an inner product on the Lie algebra. 
Given a vector field u E X(T n ) we define the corresponding momentum density by 
m = Au. The EPDiff equation for geodesic flow on Diff(T ?1 ) is given by 

(8.1) m t +£ u m = 0, 

where £ u m denotes the Lie derivative of the momentum one-form density m in the 
direction of u (see [HM] l Identifying T n ~ W n /Z n and letting x\ i = l,2,...,n 
denote standard coordinates on M. n , this can be written as 

dmi ~dmi dv? 

-W + U dxJ +m ^ + mdlv{u)= ° 

where m = midx l ®d n x and u = u l d/dx l . The CH equation is the (1 + 1) -dimensional 
version of the EPDiff equation (|8.1|) when A = 1 — A. It is natural to define a 



multidimensional generalization of uCH as equation (8.1 ) with the operator A defined 
by 



Au = (u - A)u = / ud n x - Au. 

Jrpn 
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Similarly, using higher-order tensor densities we arrive at mul tidim ensional versions 
of DP and /iDP. Let m = rriidx 1 (g> d n x <8> d n x. Then equation (8.1 ) becomes 



1.2) 



drrij 



dt dx 3 
In 1 + 1 dimensions this reduces to 



()nij 

+ u - . + rrij r + 2mjdiv(u) = 0. 



du 3 
3 dx 1 



m t + rn x u + 3mu x = 0. 



Hence, when A = 1 — A, equation (8.2) can be viewed as a multidimensional DP 



equation, while if A = /i — A it can be viewed as a multidimensional /iDP equation. 



Appendix A. Proof of Theorem 16.21 



In this appendix we prove Theorem 6.2 that is, we show that the multi-peakons 



defined in (6.8) are weak solutions in the di stribu tional sense of the ^-equation (1.2) 
if and only if \g l ,Pi}^ evolve according to (6.10). 



Equation (1.2) in weak form reada^] 
(A.l) 



1 3 — 

u t + 2 i 11 - 2 )* + x K u )a' *u-\ — - 



A 



9 * (u 2 x ) = 0, 



where g is the Green's function defined in (6.7). This is to be satisfied in the space 
of distributions V'(R x S 1 ) 



i.e. 



— oo J S 



1 



u4>t + ^u 2 4> x + \n{u){g * u) 



+ 



A 



{g * u 2 x ) 



dxdt = 0, 



for all test functions <j>{t,x) G C^(R x S 1 ). 

Let gj := g{x — q 3 ) and g'j := g'(x — q 3 ). We specify the value of g' at x = by 
setting </(0) := as in equation (6.9). We will need the following identities which 



can be verified by direct computation 

1 



(A.2) 

(A.3) 
(A.4) 



(A.5) 



9 *9j 



9j 



13 
12 



9j, 



i i 

9 *9i 



-(9i-9j)(g'i 



3 

=1 ~9j, 

--2g' J (q t )g i + 2g / i (q 3 )g j 
13 
12 



+ 9{Q i ~ Q 3 ' 



9i+ [9 



13 
12 



9i9j 



--9i + 9o + sV " <l j )(9i ~ 9') + 9(q l ~ Q 3 ) ~ 3. 



These identities hold pointwise for all x,q l ,q 3 £ S 1 (except (A.5) which holds for all 
x -,q % -,q 3 £ S 1 unless x = q l = q 3 ; a fact that does not matter in what follows). Let 
u = YliPi9i be a peakon. As explained in the appendix of |Lu] . we can compute the 



left-hand side of (A.l) as a distribution in the variable x without having to involve 



test functions explicitly. This yields 

N 
i=l 



N 



(u 2 ) x 



^2 2 PiPj9i9j- 



14 Here * denotes convolution: for two functions /, g : S 1 ~ [0, 1) 

fa f( x - y)g(y)dy- 



have (/ * g)(x) = 
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Using (A. 2) we find 



\/i(u)g' * u 



f ( 13\ , 

,7=1 V 7 



Moreover, 



X] PiPj9i9j, 



which in view of (|A.2|), (|A.3|), and (|A.5|) implies that 
t/ * (t£) = 



1 N ( 13 \ 

o X] ( W + ~ To ^ + & + ~ -* ))• 



Using these ingredients we can write ( A.l ) as 

N N 



N N A N ( 13\ 

X] - P*9*fl£) + 5^ /'-/'; - 3 Pip i [9j ~i2j 9 'j 

i=l ij=l j,j'=l 



3 - A 
6~~ 



A' 



*J=i 

We rewrite this as 



13 

7i0i + <?;<?; - ^(9i + 9j) + 3<?V - Q J )(9i ~ 9j) ) = 0. 



A-3 * 



Using (A. 4) we find 



N 



13 



X) \ \Pi + IPi^PjSfjti) ]9i + Pi\ ^2 Pi + ^2Pj9(<f ~ q j ) ~ 9* ] 5i 



i=l 



iV 



A — 3 \ ^ , 

The last sum on the left-hand side can be written as 

N 

( A ~ 3 ) PiPid'iq 1 - q J )gi, 
«j=i 



so that, using the definition (6.11) of {u x (q 1 )} and the fact that g(0) = 13/12, we 
arrive at 

AT 

£ ((Pi + (A - l)MM<f)})^ + - <f)<?-) = 0. 

i=l 

Since ^ and g[ form a linearly independent set, this equation holds if and only if 
{q i ,Pi}i evolve according to (6.10). This completes the proof of Theorem 



6.2 
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Appendix B. Proof of Theorem 16.31 



In this appendix we prove Theorem 6.3, that is, we show that the shock-peakons 



defined in (6.12) are weak solutions in the distributional sense of fiDP if and only if 



{q*, Pi, Si} i evolve according to (6.13). The objective of the proof is similar to the 
corresponding proof for DP presented in the appendix of |Lu] . However, since we 
consider the spatially periodic case and the Green's functions are very different, the 
details of the two proofs are quite different. 
The fiDP equation in weak form reads 



[B.l) 



u t + ^{u 2 ) x + 3/i(n)g' * u = 0, 



which is to be satisfied in distributional sense. We use the same notation as in the 



to verify that (B.l 



3fj,(u)g' * u are rep 



proof of Theorem 6.2 Let 5i(x) = 5(x — q 1 ) and note that g" = 1 — <5j. It is enough 
holds as a distribution in x when the functions ut, ^(u 2 ) x , and 



aced by 



and 



N 
i=l 



1 N 

« (« 2 )« = (PiPjdid'j + 8iSjg'i(l - Sj) + PiSj(g'ig'j + 9i0- ~ s j))) 



(B.2) 



3p,(u)g' * u = ^ (~PiPj {^jj - ^ g'j + 3piSj(l - gj)^j 



respectively. In order to derive (B.2) we used the fact that 



A' 



g' *u = J2 i p ^ 9 ' * + s ^ 9 ' * 9 'r>) 
i=i 



together with the identities (A. 2) and (A.3). 



Putting these expressions together we find that equation (B.l ) can be written as 



N 



^ (Pi9i ~ Pitfg'i + Si9i ~ Sitfi 1 - 



4=1 



N 



+ ^ {PiPjSiQj + SiSjQiO- - ^j) +PiSj(g'ig'j +gi( 1 - s j))) 

+ E (-Pi p 3 (dj ~ g'j + 3p4Sj(l - gj) \ = 0. 
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Employing (A.4) we can rewrite this as 



13 



E i \Pi + 2 Pi ^Pjg'jiq*) ) 9i +Pi ( + ^2,Po9{<t - q 3 ) - o 1 J g'% + kg'i - s« 



i=l 1 y j^i 
(B.3) 

N 



E {s i s j g' i +p i s j (g' i g' j + gi)) 

JV N I 

3 E ftflj(l - g s ) + E Si lq l - E s jg'j ~ ^Pjgj ] S i = 0- 

i,j=l i=l \ j j 



The identity (A. 5 ) shows that the two terms Yli jPi s j{g'ig'j~^gi) an d 3 ^ ^ piSj(l—gj) 
combine to give 



N 



E PiSj (2( 9l - 9j ) + g'(q l - - fl}) + g(q* - q>)) . 



We find that (B.3) can be written as 
E [pi + 2p i ^2 Pj g' j (q l ) + 2 Pi E Sj - 2 Si E^j I gi 



3 



N 



+ E P M ^2Pj9j((f) + E s j9 , (Q t ~ I 3 ) ~ <f +Si + Si^2 s j ~ Si^Pjg'W 



q 1 ) g'i 



*=1 \ \ 3 

N 



N 



+ 12 Si \ E^tf ~ ^) ~ ?* I + E a * K ~ E " E^* \Si = 0. 

%=i \ j J i=i \ j j 

The term Sj ( ^ ■ P jg{q l — q>) — q l \ equals 



Using that 



3 3 

K(^)> = Ew 5 '/' 1 • E- s > = E^> 

we arrive at the equation 

N 

E (Pi + 2p i {'U ;c (g i )} - 2si{u xx (q 1 )}) 5i + E (** ( n (^) - 9*) + k + a i {ti a (g < )}) ff- 



i=l 



i=l 



+ e - <n + E s * («* - * = °- 

i=l i=l 
Since {1, gi,g'i,5i} form a linearly independent set this equation holds if and only if 
{q l ,Pi, Si\i evolve according to (6.13). This completes the proof of Theorem 6.3 
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Appendix C. Conservation of H 2 for ^PjP 



In this appendix we verify explicitly that the functional H2 defined in (4.1) is 
conserved under the flow of fiDP. The following two forms of fiDP are used: 

u t + uu x + 3fi(u)A~ 2 u x = 

and its derivative 

u tx + uu xx + u 2 + 3(/i(-u)) 2 - 3ufi(u) = 0. 
We also use the identities fi(A~ 2 u) = fi(u) and A~ 2 u xx = —u + fi(u). We compute 

^ = - (?>^u){A- 2 u x ){A-^u xt ) + l -u 2 u t 



Si 



^3/i(u)(A A1 2 u x )A /t 2 (uu xx + u 2 x + 3/x(u) 2 - 3u/x(u)) 



+ ^u 2 (uu x + 3n{u)A^ 2 u x ) 

We use uu xx + u 2 = \d 2 {u 2 ) and find 
dH 2 



dt 



Si 



{-^{u) 2 {A; 2 u x )A- 2 u + MufA; 2 u x + ^{u){A- 2 d 2 x {u 2 ))A- 2 u x 



+ ^u 3 u x + ^u 2 ^(u)A^ 2 u x 



Here the first two terms and the fourth term vanish by periodicity. The third term 
is equal to 

^(u)(A- 2 d 2 (u 2 ))A- 2 u x = ^(u)(-u 2 + fi(u 2 ))A~ 2 u x 
whose first part cancels the fifth term and the second part vanishes by periodicity. 
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